Dephasing in ferromagnetic nanowires: The role of spin waves 
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We present a calculation of the dephasing time of electrons in a ferromagnet relevant for the con- 
ductance fluctuations. We focus on the contribution from the interaction with spin waves. Explicit 
results are presented for a quasi-one-dimensional systems. Going beyond previous calculations, we 
do not restrict ourselves to the limit of a small exchange splitting compared to the electronic elastic 
scattering time, nor does our calculation rely on the diffusion approximation to describe electronic 
transport. 



I. INTRODUCTION 

Effects of quantum coherence on transport proper- 
ties of disordered metals have been intensively inves- 
tigated in the past few decades, and are by now well 
understood. 'i'Sl Examples of interference phenomena are 
weak localizatiorPHSl and Aharonov-Bohm oscillations.'^ 
Another signature of quantum coherent transport is the 
phenomenon of 'universal conductance fluctuations' 
The conductance of a disordered metal exhibits random 
but reproducible fluctuations as a function of an external 
control parameter, such as the magnetic field or a gate 
voltage. The typical amplitude of these fluctuations is of 
order e'^/h at zero temperature, independent of the mi- 
croscopic details of the disorder. At finite temperature, 
the conductance fluctuations are suppressed because of 
thermal averaging and dephasing. 

Quantum coherent transport in ferromagnetic dis- 
ordered metals has attracted attention only at a 
later stage, stim ulated by the discovery of the giant 
magnetoresistanc^iSil! and by the emergence of the field 
of spintronics.'iS! In spite of the presence of an internal 
magnetic field, weak localization is observed in ferro- 
magnetic conductors if the sample size or phase-breaking 
length is small enough that the effect of the orbital mag- 
netic field or sp in-orb it scattering can be neglected on 
that length scalcli^*^ Several groups have measured con- 
ductance fluctuations in ferromagnetJ^^'^ and obtained 
estimates for the dephasing rates and their temperature 
dependence from these. '^^** In some of these experi- 
ments, the measured dephasing rates were significan tly 
larger than in otherwise comparable normal metalsj ^^"^^ 
but the microscopic mechanism responsible for this en- 
hanced dephasing could not be identified. 

With a static and spatially uniform magnetization, the 
propagation of majority electrons and minority electrons 
in the ferromagnetic metal is decoupled, and the quan- 
tum corrections to transport are essentially equal to those 
in a normal metal. Differences only appear in the pres- 
ence of spin-orbit coupling, or when the magnetization 
is not constant as a function of time or position, as is 
the case with domain walls or in the presence of spin 
waves (magnons). The weak localization correction to 
the conductivity of ferromagnets has been theoretically 



investigated in the presence of spin-orbit interaction,'^^! 
as well as in the presence of domain walls .'^21 The effect of 
domain walls or spin-orbit couplin g o n the conductance 
fluctuations has also been explored.'^^HSSI 

Spin waves, fluctuations of the magnetization direction 
that vary in time and space, are expected to contribute 
to the dephasing of electrons in a ferromagnet. Takane 
has evaluated the dephasing rate relevant for the con- 
ductance fluctuations in quasi-one-dimensional wires.'^ 
Muttalib and Wolfle performed a similar calculation for 
the phase relaxation rate in thin ferromagnetic films.'^ 
The calculation of Ref . 2Ij is restricted to the limit where 
the exchange splitting A is much smaller than the elec- 
tronic elastic scattering rate h/rd, which is typically not 
the case in (elemental) ferromagnets. Reference [2H] ad- 
dresses both the clean and diffusive limits, but still relies 
on the inequality qtIbI Ij where qt is the wave num- 
ber of a thermal spin wave, and lei = vfTbI the elastic 
mean free path. This limits the applicability of the the- 
ory to rather low temperatures, where dephasing from 
electron-electron interactions is likely to dominate over 
spin-wave-induced dephasing. 

In this work we present a calculation of the spin- wave- 
induced dephasing rate of electrons in a ferromag- 
net. We derive a general expression for I/t^ in terms of 
the spin wave dispersion relation for a model of Gaussian- 
white-noise disorder, and calculate the explicit temper- 
ature dependence for the special case of a wire geome- 
try. In contrast to previous calculations, we do not re- 
strict ourselves to the limit of small ATd/h or employ 
the diffusion approximation. We not only find qualita- 
tively different results for the clean limit ATc\/h ^ 1 
and dirty limit Atci/Ji <^ 1, but inside the clean limit we 
also find qualitatively different dephasing rates for the 
cases when majority and minority electrons have equal 
or different elastic scattering rates. In the most real- 
istic regime Arai/h 3> 1, we find that the temperature- 
dependent dephasing rate obeys a power law: at low tem- 
peratures 1/t0 oc T^/^, whereas at higher temperatures 
cx These temperature dependences are such 

that dephasing from electron-electron interactions, which 
scales cx T^/^ in a wire geometry,** is stronger than de- 
phasing from spin waves at the lowest temperatures, but 
they leave open the possibility that spin waves are the 
dominant source of dephasing at higher temperatures. 
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The structure of the paper is as follows. In Sec. |TT] 
we first present all ingredients of the model we use to 
describe the propagation of s-band electrons in the dis- 
ordered ferromagnet and their interaction with the local- 



ized spins of c?-band electrons. In Sec. Ill we describe a 
calculation of the dephasing rate in a semi-classical pic- 
ture, considering the coherent propagation of an electron 
in a fluctuating classical exchange field. In Sec. |IV| we 
then outline our full diagrammatic calculation. Details 
are kept for the Appendix. In Sec. |V] we finally calcu- 
late the temperature dependence of the dephasing rate, 
focusing for simplicity on a quasi-one-dimensional wire 
geometry. We conclude with a discussion of our results 
in Sec. 



II. MODEL 

For the conduction electrons in the disordered ferro- 
magnet, we consider the effective single-particle Hamil- 
tonian 



V- 



(1) 



where i?kin = p'^ /2m is the kinetic energy, V the impurity 
potential, and Hsd describes the exchange between the 
conduction electrons and the d electron spins. For Hsd 
we take the simple model 



Hsd = -J's(r) • <T, 



(2) 



where J is the exchange constant and s(r) is the spin 
density of the d-electrons expressed in units of h. The 
exchange constant therefore has dimensions energy times 
volume. In a ferromagnet, the magnetization s(r) fluc- 
tuates around a nonzero mean value s = se^, which we 
take along the z axis. To linear order in the fluctuations, 
only transverse fluctuations of s need to be considered, 
so that we can write Hsd as 



Hs 



1 



(3) 



where A = 2Js is the exchange splitting between major- 
ity and minority electrons and 



Hsd,±iT^,t) - J 



s_(r,0 
^{r,t) 



(4) 



where s± — Sx isy. (The exchange splitting A may be 
different from Js if exchange interactions between con- 
duction electrons are taken into account. Here and be- 
low, we will not make use of the relation A = 2 Js implied 
by our model, but consider A an independent parameter 
instead.) 

In our calculations, we will take A small in compari- 
son to the Fermi energy, so that majority and minority 
electrons have the same density of states f at the Fermi 
level and the same Fermi velocity vf- In order to model 
the effect that majority and minority electrons have dif- 
ferent scattering rates in realistic ferromagnets, we take 



different disorder potentials for majority and minority 
electrons, 



V = 



_ I 



V, 



(5) 



Here Vf and Vi are Gaussian white noise potentials with 
correlation functions 
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5(r-r'), 



where and are the elastic mean free times for major- 
ity and minority electrons, respectively, and Td is a time 
that describes the degree of correlation between the effec- 
tive impurity potentials for majority and minority elec- 
trons. The case of equal impurity potentials for majority 
and minority electrons corresponds to the case — 
and Td — ^ oo. 

The alternative to our implementation of different scat- 
tering rates for majority and minority electrons through 
a spin-dependent impurity potential is to take the Fermi 
energy difference between majority electrons and minor- 
ity electrons seriously. However, taking A to be compara- 
ble to the Fermi energy ep is incompatible with the semi- 
classical limiting procedure at the foundation of the dia- 
grammatic calculation of the dephasing time,!^ because it 
no longer allows one to separate diagrams that are small 
by a factor Ti/e-pT (which are usually neglected) from 
those that are small by a factor Ti/I^t (which need to 
be kept for a diagrammatic calculation of the dephasing 
rate) . The choice of a spin-dependent impurity potential 
circumvents these problems.^^^ 

In second-quantized language, the transverse sd ex- 
change Hamiltonian Hsd,± takes the form 



Hs. 



(7) 

with Sq -t — J drs±{r,t)e "i "" the Fourier transform of 
the spin density and V the volume of the ferromagnet. 
Dynamical processes involving the excitation and absorp- 
tion of a d-band spin wave are characterized by the sus- 
ceptibility 



X^(q,r) 



ie(r)([sq,_(r),Sq.+ (0)]), (8) 



where 0(t) = 1 for r > and 0(t) otherwise is the 
Heaviside step function. The susceptibility (Qi ''') de- 
scribes the response of the d-electron spin density to an 
applied magnetic field. The Fourier transform X-+ (Qi ^) 
is conveniently expressed in terms of the spin wave fre- 
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quencies uj^^MM 



2s 



(9) 
(10) 



Lu — LU^ + ir] 
where 77 is a positive infinitesimal. 

III. SEMICLASSICAL PICTURE 

The spin wave contribution to the dephasing rate in 
ferromagnetic metals can be understood from a semiclas- 
sical argument. Hereto we consider the d-electron spin 
density s(r) as a classical variable, and look at the ef- 
fect of fluctuations of s(r) on the coherent propagation 
of conduction electrons. We treat the orbital degrees 
of freedom classically, with conduction electrons moving 
along classical trajectories r{t). Their time evolution is 
governed by the sd-exchange Hamiltonian ([2|, which is 
taken to act on the spinor degrees of freedom of the con- 
duction electrons only. In this approach, the Hamilto- 
nian Hsd depends on time t explicitly through the time 
dependence of the d-electron spin density s and implicitly 
through the time dependence of the conduction electron's 
position r(t), combined with the r dependence of s. 

As in the previous section, we separate Hgd into a 
(large) time-independent part Hsd = ~(A/2)crz and a 
(small) time-dependent part Hsd,± describing the trans- 
verse fluctuations of the d-electron spin density. We 
describe the time-evolution of a spin described by the 
Hamiltonian Hgd in the interaction picture, in which 
Hsd,± is treated as a perturbation. The time-evolution 
operator reads 



(11) 



where the perturbation is written in the interaction pic- 
ture. 



to consider the diagonal elements U^^{t) and f/^^(r). We 
also note, that in the expansion of the the diagonal ele- 
ment i/j^^ (r; t) the transverse spin densities always occur 
pairwise in the combination s_(r[t2], i2)s+(i"[ii], ii), with 
t2 > ti. For T much longer than the time scale for fluc- 
tuations of the magnetization, these events appear well 
separated in time, and one finds that the leading contri- 
bution to C//_-|-|-(r) is given by 



(14) 



where 



i^_+(r2,Tl) = 3^{r[t + T2],t + T2)s+{r[t + Tl],t + Tl) 

X e(T2 - ri)e-*^("^-"i)/'*. (15) 

The averaging bar denotes an average with respect 
to t. For long times, the magnitude of J7/,^^(t) decays 
cx exp(— r/2T0_^), where T^^-f is the dephasing time. From 
Eq. ( 14 1 we then conclude that the dephasing time r^^^ 
is given by 



1 



t2 roo 

Re / dT'{F_+iTi+T',Ti))r„ 

n Jo 



(16) 



where the brackets (•••) indicate an average along the 
trajectory r{t + ti). Repeating the same analysis for 
minority electrons gives 



1 



2 Re / dr'(F+_(Ti-Hr',ri 
n Jo 



(17) 



with 



i^+_(r2,Tl) = S+{r[t + T2],t + T2)s^{r[t + Tl],t + Tl) 

X e(T2 - ri)e*^("=-"i)/'*. (18) 



It remains to evaluate the time averages in Eqs. ( 15 ) 



and ( 18 ) and to perform the average along the trajectory 
in Eqs. (16) and ( [l7[ ). The spin correlation function re- 
quired for the calculation of F ^. follows directly from 



^ldV(^) = e'-^^"^/'^^Hsd.At)^~'^^"'^'^^^ (12) the fluctuation-dissipation theorem,l2ni 



and T denotes the time-ordering operator. The rela- 
tion between J7^^^ and the time-evolution operator in the 
Schrodinger picture U is 



U{T-t) = e'(*+-)A^./2nf;(/)(^.^)g-»tA../ 



2h 



(13) 



We now calculate the average U (r) with respect to the 
time-dependent fluctuations of the transverse spin den- 
sity s± . In order to separate the explicit time dependence 
of s±{v[t],t) through the time argument t and the im- 
plicit time dependence through the t dependence of the 
position r, we consider a fixed trajectory, but at different 
starting times t. We note that the off-diagonal elements 
of the evolution operator vanish upon taking the time- 
average, because they contain an odd power of the trans- 
verse d-electron spin densities s±, so that it is sufficient 



xlm{x-+(t^,q)}, (19) 

where x h is the magnetic susceptibility, which describes 

the magnetization response to an applied magnetic field 
(note that we have set for convenience fc^ = 1). The 
imaginary part of the response function x hi'^i q) is de- 



termined by the spin wave spectrum, see Eq. ( 10 ) 



lm{x-+{u,q)} = 2tts6{lo - Lol^. 



(20) 



For the average along the trajectory r{t) we consider 
s-wave impurity scattering with elastic mean free time 
''oi = T-f = Ti and ballistic propagation between scatter- 
ing events. In Fig.llk we illustrated a part of a trajectory 
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p'{r) + / P{j)J'-^ +... 

x;<--'' X---'"" 

n = n = 1 

FIG. 1: (Color online) Illustration of our semiclassical calcu- 
lation, (a) A spin-up electron propagates ballistically through 
the ferromagnet, its direction of propagation regularly being 
randomized by scattering off impurities (black circles). At 
time t, it interacts with the fluctuating field of the d-band 
electrons and its spin is flipped down. After a time r, the 
electron has traveled over a distance p and it interacts again 
with the d-band electrons, its spin being flipped up again. 
Since the interactions occur randomly in time and space, the 
phase of the electron becomes randomized while propagat- 
ing, (b) The part of the trajectory important for dephasing 
is when the electron carries spin down. When averaging over 
the displacement p(t) during this part, we sum over trajec- 
tories with no scattering events (n = 0), with one scattering 
event (n = 1), etc. 



For small q and Tiuj^^ <Si A, this expression takes the form 
f dq AJ^Ts Dq^ 

J (2^)3 hu-^ A2[l + (ATel/;i)2]' ^ ^ 

where D — WpTei/3 is the diffusion constant. 



IV. DIAGRAMMATIC CALCULATION OF 
DEPHASING TIME 

We now outline a full diagrammatic calculation of the 
spin wave contribution to the dephasing time. In the 
language of diagrammatic perturbation theory, the prop- 
agation of electrons in a disordered metal is described by 
means of the diffuson and Cooperon propagators. These 
describe the phase-coherent propagation of multiply scat- 
tered electrons traveling along identical paths or along 
identical but time-reversed paths, respectively. In the 
case of a ferromagnet, the internal magnetic field re- 
sults in rapid dephasing of electrons traveling along time- 
reversed paths, and the Cooperon is strongly suppressed. 
For this reason, the calculation below focuses on the dif- 
fuson propagator. 

The diffuson propagator relevant for the calculation of 
conductance fluctuations reads 



for a spin-up electron: The electron scatters off impuri- 
ties (black circles), and at times t and i -I- t it interacts 
with the fluctuating field of the d-band spins causing a 
spin flip. The electron thus carries spin down for a time 
r (the red dotted part of the trajectory), and during this 
time it travels over a distance p. 

The average over trajectories in Eqs. (16 1 and (17) 

iq-p(r) 



amounts to averaging 



over all possible dis- 



placements p in time r. Summing over all possible num- 
ber n of scattering events in the time interval r (see Fig. 
[T]d), we obtain 



1 

T0 



Re 



OO 



dq U'^Ts 



q 



dq AJ^Ts 



Re 



2qvpTei — iK 



(21) 



where we have dropped the spin index for t^, because 
minority electrons and majority electrons have equal de- 
phasing times in the semiclassical model, and 



K = -2tqvF / dre-'^^/'^+'^D^-^/^ei 
Jo 

dO smg -^^^^^p^g 



In 



2 
1 + l{^/fl + ~ qVY)T,i 
1 + i(A/?i + u;^« + 9UF)Tcr 



(22) 



(24) 



r', r, £ 



huj/2)G^,{v,v',e-hLo/2)), 



where the bar denotes an average with respect to 
quantum-mechanical and thermal fluctuations of the ri- 
band electron spin density, whereas the brackets (• ■ ■ ) de- 
note a disorder average with respect to the impurity po- 
tential V . Further and are retarded and advanced 
Green functions for the conduction electrons, respec- 
tively. When calculating conductance fluctuations, the 
two Green functions in (25) are part of interfering elec- 
tronic trajectories at significantly different times. Since 
interactions with the spin density fluctuations along the 
two different trajectories are thus completely uncorre- 
lated, the averaging over the fluctuations is done sepa- 
rately for the two Green functions. 

We are interested in the Fourier transform 



I?,,,<(q;e,u;) = J dpV,^,,{p;e,uj)e'''f . 



(25) 



Since the dominant contribution to coherent propagation 
comes from constructively interfering trajectories of elec- 
trons with the same spin (due to the exchange splitting, 
electrons with different spin dephase rapidly) , we will fo- 
cus on the diagonal elements I?(j^(j(q; e, w). For small q 
and uj this diagonal propagator has the asymptotic pa- 
rameter dependence 



2?<j„(q;£,a;) 



2'KV 



1 



(26) 
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where = VpTa-/3 is the diffusion constant. This pa- 
rameter dependence corresponds to an exponential decay 
a exp(— r/r^^cr) of the Fourier transform of V^^cri^.] £, ^) 
at q = 0, similar to the exponential decay of the ab- 
solute value of the evolution matrix element Uaa in the 
semiclassical picture of the previous Section. 

We denote the diffuson propagator in the absence of 
coupling to transverse fluctuations of the d electron spin 

density by ^^^^(q; e, w). This "bare" propagator is given 
by the equation 

2'^i(q;e,'^) = uyr, (27) 



l-ni°i(q;e,c^) 



where 



n(°l(q;£,c.) 



(28) 



5](G^(k, s + hu/2)) (G^(k - q, £ - hLu/2)) 



is the bare structure factor, see Fig. [2]i. Substituting 

1 



(G?(k,£)) 
(G^(k,£)) 



e - £fc + M + Aa/2 + hi/2T^ ' 
1 

e - Efc + M + - h.i/2T„ ' 



(29) 
(30) 



for the impurity-averaged single-particle retarded and ad- 
vanced Green functions, one finds the diagonal element 



In 



(31) 



in the limit q <^ fcp, h\uj\ <^ fi, which leads to the asymp- 
totic small-g and small-w dependence of Eq. ( 26 1 without 
the dephasing term. 

Coupling to spin waves alters the structure factor, ulti- 
mately leading to a finite dephasing time in the diffuson 
propagator (251. A diagrammatic evaluation, details of 



which can be found in the Appendix, leads to two types 
of corrections to the structure factor. The first correction 
is a spin-wave-mediated renormalization of the diffusion 
constant, which does not alter the functional form of the 
low-q and low-w asymptotics of the diffusion propagator. 
This correction is the equivalent of the Altshuler-Aronov 
correction to the conductivity from electron-electron in- 
teractions. Since this is a completely elastic correction, 
which does not contribute to dephasing, it will not be 
discussed here. The second correction is the spin- wave- 
mediated contribution to the dephasing rate. The rele- 
vant diagrams for this correction are shown in Fig. [2J3. 
Denoting this contribution by Il™{q;e,uj) and analyzing 
the asymptotic dependence of the diffuson propagator for 
small q and small lu, one finds that 



(0;e,0) 



(32) 



(a) e+^hij k a 
: 



£—^hui k— q a 



(b) 



k cr r l . k cr 



k-Qs, 



sw 







~> r 

:_<_^_L 



e — \ huj k — q a 



X XX 

= -< — — 



FIG. 2; (Color online) (a) Diagrammatic representation of 
the bare structure factor 11^'^^. The solid blue (dashed red) 
line represents a retarded (advanced) Green function, and the 
dotted line with the cross a correlated impurity scattering, (b) 
Diagrams representing the correction to the structure factor 
from spin wave induced dephasing. Correlated excitation and 
absorption of a spin wave changes temporarily the energy, 
momentum and spin of one of the propagators. The shaded 
block represents a ladder of impurity scatterings. 



Evaluation of the diagrams in Fig. [2]d yields 
,0) 



n™(0;e,0) 
_ 2J2 



X coth — — — h tanh — — 

\ 2T 2T 

X Re |n^°^^(qsw;e- ?i'^sw/2, Wsw, 2) 

^ rCTn^°'_^(qsw;£ - ?iwsw/2,a;sw, 1)^ 
To - r<:,n^°^^(qsw;e - ?ia;sw/2,a;sw,0) 

where we abbreviated 

ni°-^(q;e,w;n) 



,(33) 



27rt/T^t/ 



^(G^(k,e + ?.c./2)) 



X (G^(k, e + hoj/2)r{G^,{-k q,e- nLo/2)) 



2h"'qvF 1 — i{uj + Aa /h + qv-p)Ts^ 



with 



1 1 



1 



1 



To 2t^ 2t^ Td' 



2Tt 2t; 



(34) 



(35) 
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the scattering rate for correlated scattering events of ma- 
jority and minority electrons and the mean scattering 
rate, respectively. Setting £ = 0, we thus find that the 
spin wave contribution to the dephasing rate of conduc- 
tion electrons at the Fermi level is 



previous calcul ation s of the spin wave contribution to the 
dephasing ratej^SEB the expression derived above neither 
relies on the diffusion approximation ql^i ^ 1, where /ei 
is the elastic mean free path of the conduction electrons, 
nor on the dirty-limit condition Ard/h <^ 1. 



1 



1 



ToiK 



^ Re , 

V sinhihu^J' /T) 2qvFTo - iK 



(36) 



V. TEMPERATURE DEPENDENCE OF 
DEPHASING TIME IN A WIRE GEOMETRY 



with 



if = In 



i + i{A/n- 



qvF)Ta 



(37) 



This result is identical to the result of the semiclassical 
calculation, up to the replacement of the summation over 
q by an integral and the replacements sinh(?ici;q^/T) — >■ 
/T and ro,ra — Td. We note that, in contrast to 



The only ingredient still missing is the spin wave dis- 
persion relation ui™. As long as the typical wave length 
of the spin waves is large compared to the Fermi wave 
length, macroscopic spin wave theory can be employed 
to find uj™. For definiteness, we assume the anisotropy 
to be uniaxial and restrict ourselves to a quasi-one- 
dimensional (wire) geometry. Choosing the z-axis to be 
parallel with the easy axis of the magnet, one then ha^i' 



2gnBK 



AnjMs sin^ 6^ 



2gnBK 



(38) 



Here, D™ is the spin wave stiffness (usually of the order 
]J^^ A/hkp), gfiBB is the electronic Zeeman split- 
ting due to an externally applied magnetic field along 
the z axis, AI^ — gfJ-ss is the saturation magnetization, 
7 = g^BlJ'Q is the gyromagnetic ratio, K is the energy 
density characterizing the magnetocrystalline anisotropy, 
fiQ is the permeability of free space, and denotes the 
angle between q and the z-axis. In general, the demag- 
netizing field should be added to the applied field, cor- 
responding to the replacement B ^ B — /ioAfs^/47r in 
the equation above, where ^ is a dimensionless constant 
determined by the shape of the ferromagnet, such that 
^ is of order unity for a reasonably symmetric shape, 
and ^ — >■ in the limit of a quasi-one-dimensional ge- 
ometry (length much longer than width W). [Note that 
this dispersion relation differs from the one used in Ref. 
UHl ?iWq™ = Eg + Aq^. This (more phenomenological) 
dispersion relation has only two ingredients, a spin wave 
gap Eg and a stiffness A, covering only a subset of the 
regimes we investigate in this work.] 

Below, we will neglect the contribution of the 
anisotropy to the spin wave energy, supporting this as- 
sumption with typical parameters for iron {K « 5 x 
10* J/m3, 7M, « 0.25 meV, see Ref. giving 
K/2Tr^QM^ « 2 • 10^^. Then, for the case of no exter- 



Spin waves with an energy huj'^ > T do not signifi- 



nally applied field, the dispersion relation (38) simplifies 
to 



hD'^^q^ 



47r7AfsSin'^ 0q)hD^'"q'^. 



(39) 



cantly contribute to dephasing. From Eq. ( 39 ) we then 



conclude that the order of magnitude of the largest mo- 
menta which have to be taken into account in the calcu- 
lation of the dephasing rate satisfies 



(40) 



The calculation of the previous section is thus valid as 
long as (/max ^ kp, which is satisfied for all T ^ A. 
A stronger constraint has to be met in order to justify 
the use of the diffusion approximation for the conduction 
electrons, as employed in Refs[?7land[251 This constraint 
is Qmax ^ Ici^j which is met only if T <C A(A:f^oi)~^, im- 
posing a severe restriction on the range of temperatures 
which can be investigated. 



Eqs ( 36 )-( 38 ) form the most general result of our work: 
They allow to calculate (numerically) the dephasing time 
for arbitrary A (as long as A ^ T) and, by including 
a demagnetizing field into ( 38 ) , for arbitrarily shaped 



samples. Let us now investigate different limits of Ard/h, 
in which the expression for the spin wave contribution to 
the dephasing rate simplifies considerably and we can 
arrive at explicit expressions for the dephasing rate. 



A. The limit Ard/Ti ^ 1, and 7^ 

We first consider the most realistic case, i.e., where 
the exchange splitting A is large enough that AT^/h ^ 
1 and where the elastic scattering times and are 
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significantly different for the two spin directions. In this 
regime, we expand Eq. (36) to leading order in ?i/Atci, 



1 



1 



T^,, V ^ AVd sinh(?iw--/T) ■ 



(41) 



In the calculations that follow, we will replace the sum- 



mation in Eq. (41) by an integral and use the approxi- 
mation (39) for the spin wave dispersion relation. This 



replacement and the use of this approximation requires 
a further discussion for two reasons, (i) The sum in Eq. 
( 41 ) diverges for small in a wire geometry if we use Eq. 



( 39 ) for w™. The origin of the divergence is that — > 
for g — >■ at zero magnetic field when the anisotropy K 
and demagnetizing factor ^ are set to zero. The diver- 
gence is removed if we consider the full dispersion re- 



lation of Eq. (38), for which takes the finite value 



= ijJk — zg/iB^/?i-Ms at g = (and zero magnetic 
field). For sufficiently large anisotropy ujk ^ D^^/W'^, 
the integrand has a smooth dependence on q and the sum 
over q may be replaced by an integral. The resulting in- 
tegral, however, turns out to be largely independent of 
K as long as <C T, which justifies the use of a spin 
wave dispersion relation with K = in combination with 



the replacement of the summation in Eq. (41 ) by an inte- 
gral in the parameter regime hD™ /W^ < Tllok ^ T. (ii) 
The replacement of the summation over q by an integral 
is allowed only if the transverse dimensions of the wire 
W are small enough that W » q^l^ ^ kp'^^/KjT. Tak- 
ing realistic parameters for Fe and a wire width W = 20 
nm, this limits the following calculations to temperatures 
T > 0.08 K. 

We now focus on different limits where we can explic- 
itly solve this integral. The result of the integral still 
depends on the ratio between the temperature and the 
demagnetization energy AnjMs. When we express the 
demagnetization energy in terms of temperature (e.g., 
iiT'jMs ^ 3 meV ~ 35 K for iron), we see that both 
the low-temperature limit T ^ Air-fMs and the high- 
temperature limit T 3> AttjMs could be experimentally 
relevant. 



We start by rewriting the integral version of ( 41 ) as 



1 A r^\n ^'sinfl 

^ A dx dO ^ 

''''!>,<' Jo Jo sinh (ax\/x'^+ sin^ 6 



(42) 



with A = (2j2s/7r2AVd)(47r7M^/ni:)"*)3/2 and a = 
Att^Ms/T . In the low-temperature limit we have a 3> 1, 
and we can approximate sin 6 6 for all relevant 6. We 
then can perform the integration, to arrive at 

1 /2t5/2 

with Ci « 0.0473 being a numerical prefactor. The op- 
posite case of large temperature has a <C 1, and can 
be evaluated as well. In this case we can approximate 



■ sin^ « a; in the integrand of ( 42 ) . The resulting 



integral can be solved, yielding 

T^.. " 'A2Td(?ii^«-)3/2' 



(44) 



where the numerical prefactor C2 ~ 0.303. 

Another limit which can be considered is that of a 
large external magnetic field (but still A ^ qix-qB). 
When qiibB ^ An^Mg, we can approximate huj^ k, 
TiD^^q^ + qii-qB^ independent of 6'q. For high tempera- 
tures T 3> giJ-sB we then recover the high-temperature 
result of Eq. (44). In the low-temperature limit T <C 



gfisB all spin wave energies contributing to the inte- 
gral in Eq. (36) are larger than T, so we approximate 



l/sinh(?iw^'^/r) w 26"'*"^/^. The resulting integral can 
be calculated, yielding a dephasing rate that is exponen- 
tially suppressed with increasing magnetic field B, 

— = — J^sT^^^ g-ffMBS/T y o 

(45) 



B. The limit Ard/h > 1, and = 

For clean ferromagnets, for which ATd/fi is still large, 
but not large enough to cause the elastic scattering times 
for majority and minority electrons to be significantly 
different from each other, we can use = = Tq = Ta = 
Toi and Td 00. We see that in this case the leading order 
contribution we found in Eq. (41 1 vanishes and we have 
to expand Eq. (36) to higher orders in Ti/Atc\. 



Again 

limiting our calculation to the regime 

the leading-order contribution to the dephasing rate is 

found to be 



1 



dq 

(27r)3 



(Atoi)4 sinh(?ia;^"/T)' 



(46) 



where D — WpTci/3 is the conduction electron diffusion 
constant. Asymptotic expressions for the dephasing rate 
can then be obtained in the same limits as treated in 
the previous subsection. In the low-temperature limit 
T < Att^Ms, we find 

1 ri^ P DT'^/'^ 

T^.a A-*r_fj75/iB(ft£'''")^/2 

with C3 K, 0.0367. In the opposite case of large temper- 
ature T > in^Ms we find 



1 



C4 



A4T2(nL»-w)5/2' 



T:^ An J Ms, 



(48) 



where the numerical prefactor C4 « 0.297. Finally, with 
a large external magnetic field, QHbB ^ Air-fMs, and low 
temperatures, T <C g^^B, we arrive at 

1 3 n'j^DsT^/^ 



27r3/2 AV^j (?!£)'' 



')5/2 



(49) 
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We note that the three exphcit results in this Section 
are all a factor ~ {n/ATei){T/A){D/D^^) smaller than 
the corresponding results in the previous Section. The 
main reason for the large enhancement of dephasing when 
^ Ti lies in the fact that in this case electrons with 
different spins effectively see a different impurity poten- 
tial. Constructive interference of propagating electrons 
with different spin is then not only suppressed by their 
energy difference A, but also by the fact that their actual 
trajectories quickly diverge. 



C. The limit Ard/h < 1 

Finally, we consider the limit ATd/h ^ 1 of a dirty 
ferromagnet. We also assume here that A is too small 
to cause a significant difference between t-^ and r^, so 



we again set ^ = tq = Ta = Td and — )■ cx3. If 
W '3> max[kp^ A/T, h / Atc\] the summation over 
q can be replaced by an integral and we find 



1 



dq 



AJ^STelDq'^Tc 



(27r)3 [{hDq^T^i)^ + (At^i)^] smh{huj^^ /T)' 

(50) 

This integral is difficult to evaluate; The most accurate 
way to find a dephasing time in a specific material is to 
solve it numerically. We can however arrive at an order- 
of-magnitude estimate and find the power-law depen- 
dence on T in different regimes. Since Ard/h 1 and 
T <IC A is a requirement for the validity of spin wave dis- 
persion relation ( 38 ) , we only treat the low-temperature 
limit T <C Att^Ms . Performing an analysis similar to the 
one leading to Eq. ( 43 ) , we find 



I J'^D^^A^/'^ 1 



dx 2x3/2arcoth(e^) + V^[Li2(e-^) - Li2(-e-^)] 
27r3 + c2 



(51) 



with c = AD^/TD - (A/r)(ATei/n)(l/fcF/ei)^ and Lis 
the dilogarithmic function. Evaluating the integral for 
c <C 1 and c ^ I we find 



1 



Q 



75Mb A2 {HD^ 



C < 1, 



O 1, 



(52) 
(53) 



0.0514 numerical constants, 
limit ATc\/h ^ 1 has been 



with Cs ~ 0.0884 and Cg ft 
We note that the same 
considered previously by Takane.^^ Our result (50) in- 
deed coincides with Eq. (38) from Ref. 27 in the imit 
of A ^ huj™ and = (up to a factor I/tt, which 



Takane overlooked in his last step). Our estimates (52) 



and ( 53 ) differ qualitatively from those in Ref. ,27|because 
ours have been derived consistently in the limit where 
and Ti are set equal, whereas Ref. [17] considers the case 
of different scattering times for majority electrons and 
minority electrons. 



VI. CONCLUSION 

We calculated the contribution of the interaction be- 
tween s-band electrons and c?-band spin waves in fer- 
romagnets to the dephasing rate of coherent propaga- 
tion of the elec trons . Our work, as opposed to previ- 
ous calculations ,'^^'21! neither relies on the diffusion ap- 
proximation qlci <^ 1, nor on the dirty-limit condition 
Atci/Tl <C 1. We found qualitatively different results de- 
pending on whether Atci/Ti ^ 1 or Ard/Ti <C 1. 



We investigated the explicit temperature dependence 
of the dephasing rate in quasi-one-dimensional systems. 
In the most realistic limit of Atci/Ti ^ 1 we find that 
the temperature-dependent dephasing rate obeys a power 
law: at low temperatures it is proportional to T^/^, 
whereas at higher temperatures it becomes proportional 
to T"^/^. If ATc\/fi is still large, but A is not large enough 
to cause a significant difference between the effective dis- 
order potential seen by electrons with spin up and down, 
the dephasing rate becomes proportional to T^/^ at low 
temperatures and to T^l"^ at high temperatures. We also 
provide expressions for general Ard/h and arbitrarily 
shaped samples. In this general case however, the sum- 
mation over q can not be performed in closed form, and 
our general expression has to be evaluated numerically. 

A relevant question to answer is how the phase relax- 
ation time thus found compares to the phase relaxation 
time due to electron-electron interactions. For a quasi- 
one-dimensional geometry, the latter is of the order^ 



T 



2/3 



1 



T 



klAh 



2/3 



(54) 



where vi is the effective one-dimensional density of states 
and A is the cross-sectional area of the wire. Since 
dephasing from electron-electron interactions scales oc 
r2/3, which is a lower power of T than any of the low- 
temperature limits considered in Sec. |V] Coulomb inter- 
actions provide the dominant source of dephasing in the 
limit of sufficiently low temperatures. On the other hand, 
for sufficiently high temperatures, dephasing from spin 
waves may take over. 
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We now estimate the temperature Tmin above which 
the spin-wave-induced dephasing dominates for a wire 
with a cross section A = (20 nm)^, using reaUstic param- 
eters for Fe. Hereto, we take A ^ 9x 10'^ K (using param- 
eters from Ref. 33 and assuming that roughly A ~ Js) 
and Tel 3 X 10^^'' s ^ (250 K)~^, which corresponds to 
the clean limit Ard/Tl ^ 1, for which the results of Sec. 
V A apply. We use our low-temperature result ( 43 1 , esti- 
mate D^'^ ~ A/hkp and assume Td~ Tei, which, together 
with typical parameters for ironj ^^ l '^^ l gives Ty^in ^ 3.5 K. 

We have to make several consistency checks here: (i) 
The temperature found is consistent with our choice for 
the low-temperature regime T <^ Att^/Ms required for the 
validity of Eq. ( 43 ) . (ii) The transverse dimensions of the 



wire should exceed kp^ y^~AjT]^ ^ 3 nm, which is indeed 
the case, (iii) To see whether the phase relaxation time 
at this transition temperature is still compatible with 



the quasi-one-dimensional limit of Eq. (54), we evaluate 
the dephasing length ^-Dr^^ee with the numbers used at 
temperature T ~ 3.5 K: We find ^JDt^^ ~ 10 /.tm, 
indeed much larger than the chosen wire thickness of ~ 
20 nm. Although the precise value of Tmin depends on 
effects not taken into account here, such as the detailed 
band structure of Fe, or specifics of the impurities, the 
order-of-magnitude estimate Tmin '-^ 3.5 K leaves room 
for a significant parameter regime where the interaction 
with spin waves could be the dominant mechanism of 
electronic dephasing. 

Comparing our results with the available 
experiments,'^^"2II we make the following observa- 
tions: (i) The anomalously short dephasing lengths 
observed in some of the experiments (Refs I16I19I211 
report ~ 100-300 nm at 30-80 mK, and ~ 10 nm 
at temperatures around 1 K) differ from the estimated 
rate for spin-wave-induced dephasing in Fe by three 
orders of magnitude. It seems unlikely that this large 
difference can be attributed to material-specific details 
not taken into account in the theory. Hence, we believe 
that spin waves are not the source of the anomalously 
high dephasing rates seen in these experiments, (ii) The 
temperature dependence of I/t^ from experiment s on 
quasi-one-dimensional ferromagnetic metal wireJ^^'^I] 
could fit a T'^/^-dependence in some data sets. A 
problem with extracting this power law dependence from 
the data is that the value of the dephasing length is so 
short (see above), that becomes comparable to the 
transverse dimensions of the wires in a significant range 
of temperatures. This makes it difficult to convert the 
measured quantity (the rms value of the conductance 
fluctuations or the correlation field or bias voltage) 
unambiguously into a temperature dependence of L^, 
since the relation between these quantities and 
depends q ualitat ively on the effective dimensionality of 
the system .1^21251 (jij) xhe measured indicators of quan- 
tum coherence (such as noise power of time-dependent 
conductance fluctuations,'^ or correlation bias voltag^^ 
seem to be largely invariant when the external magnetic 
field is increased over a range of several Tesla. This, 



too, seems to indicate that the interaction with spin 
waves is not the dominant mechanism of dephasing 
at low temperatures, since one expects an exponential 
suppression of the spin-wave-induced dephasing rate 
when g/is-B ^ T. However, at temperatures T > Tmin, 
where we theoretically expect spin waves to start playing 
a role (we estimated Tmin ~ 3.5 K, see above, which cor- 
responds to magnetic fields ~ 2.6 T), one indeed would 
need a field of several Tesla before any field-induced 
suppression of the effect would become noticeable. 

We close with a remark comparing the present calcu- 
lation and the calculation of the dephasing time from 
electron-electron interactions.^ There are two important 
differences between these two calculations. First, for 
Coulomb interactions the dephasing rate for a quasi- 
onc-dimcnsional geometry can not be obtained in per- 
turbation theory, but instead needs a self-consistency ar- 
gument to cut off an otherwise divergent integral. No 
such divergence occurred in the present calculation. Sec- 
ond, the calculation of the dephasing rate from electron- 
electron interactions that enters the conductance fluctu- 
ations requires the calculation of a more complicated ob- 
ject than the simple diffuson propagator we consider.l^ 
The origin of the difference is twofold: (i) The Coulomb 
interaction becomes effectively long range at low frequen- 
cies w, whereas the spin- wave- mediated interaction re- 
mains short range even at small uj because of the ex- 
change splitting A. Indeed, since every interaction with 
the spin waves involves an electronic spin flip, correla- 
tions over distances larger than min[?iwF/A, (?iZ?/A)^/^] 
are strongly suppressed by the exchange splitting, (ii) 
The divergence at q = and the resulting = in 
our calculation (corresponding to a static magnetization 
pointing in a random direction) is cut off by the magne- 
tocrystalline anisotropy, physically pinning the magneti- 
zation to a fixed direction. 
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Appendix A: Evaluation of the diagrams 

The diffuson propagator describes the coherent prop- 
agation of two quantum mechanical probability ampli- 
tudes. The effect of spin waves on the propagation is 
qualitatively different for the cases that these amplitudes 
describe propagation through the sample at the same 
time, or at different times. The former case is relevant, 
e.g., for a calculation of the disorder-averaged conduc- 
tance, whereas the latter scenario is relevant, e^., for 
the calculation of the conductance fluctuations.^ Dia- 
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(a) 



< 



(b) 



FIG. 3: (Color online) Diagrams representing the spin wave 
correction to the diffusion propagator. (Impurity lines are 
not shown in the diagrams in this figure.) (a) Diagrams in 
which the spin wave propagator does not connect the retarded 
and advanced propagators in the diffuson. (b) Diagrams in 
which retarded and advanced propagators are connected by a 
spin wave propagator. Only diagrams shown in part (a) are 
relevant for the calculation of the dephasing rate. 



grammatically, the difference between the two cases is 
that all spin-wave-mediated interaction lines have to be 
included in the former case (diagrams of Fig.[3]a and b), 
whereas in the latter case only interaction lines that con- 
nect the retarded propagator to itself or that connect the 
advanced propagator to itself have to be considered (di- 
agrams of Fig. |3^) . Impurity lines, on the other hand, 
can connect all Green functions. Our calculation of the 
dephasing rate applies to the second scenario. This is 
consistent with the definition of the diffuson propagator 
in Eq. (251, where the spin wave averages are taken for 



each Green function separately and the disorder average 
is taken of the product of the two Green functions. (We 



note that current conservation protects the pole of the 
diffusion propagator in the former case.) 

For our calculation of the diffuson propagator, we first 
consider the most general quantity 



= (^CTi.<T2(i"',r,ie„)5o.3,<^Jr',r,i£„)), (Al) 

where e„ > and Sm < are fermionic Matsubara fre- 
quencies. Because the spin wave average is taken for each 
Green function separately, the Green functions them- 
selves contain no off-diagonal elements, i.e., cti = (72 and 
— o'i, so we write the diffuson as 



= {gair',r,ien)Ga'{r',r,is„i)). 



(A2) 



We further focus on the dominating contribution where 
a = a' , i.e., the diagonal elements of the diffusion propa- 
gator. These are then found to obey the Dyson equation 

2^(7,(7 (r' - r;«£n,«em) = 2TnyT^Ila,cr{r' - r;ien,i£m) 

(A3) 

where no-,cr(r' — r;i£n,iSm) is the full structure factor. 
Fourier transforming Eq. ( A3 1 , we find 



2?cr,o-(q;ien,»em) = 2nvT„Il„^cr{q;ien,iem) (A4) 
X [l + (27ri/r„)~^2?^^^(q;ie„,ze™)] . 

For the calculation of the dephasing time, we need the 
simultaneous limits q — >■ 0, e„ — £-|-iO, and — ^ e — iO. 

The diagrams contributing to n™^(q; i£„, ie„) are 
shown in Fig. |4] Evaluating these diagrams yields 
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FIG. 4: (Color online) (a) Diagrams representing the spin wave correction to the disorder-averaged structure factor, (b) 
Diagrammatic definition of the Hikami box appearing in some of the diagrams of part (a), (c) Definition of the dressed vertices 
appearing in some of the diagrams of part (a). 



J^T ^ 1 



+ 



■^0 - "^'T 27rt/r^y Sk^-'^C^ ~ qsw,i£n " (k - q, IE™ ) 



H — ^ V — V'x-+(qsw,«i's)i — V't;CT(k,«e„)t;^(k-q,ie„)2t;_^(k-q-qsw,jem-«t^.s) (A5) 



'^A 2-nvr„v a^-,ien)Q - q- q^w, - iv s)Q cr{}^ - q,ie™)]' 



To 



J^T 1 

Vs<^ryi qsw 

J2T ^ 1 



To 



^£rj 



To 



To - To 



27riyT„V' 



The impurity averaged Green functions Q^i^^ie) and response function x hC^j^^s) ^ written in Matsubara 

representation. The contribution to ( A5 ) by the Hikami boxes (see Fig. |3j)) is 



\^ — '^—Tr ^cr(k, ie„)t/o-(k - q, i£„i)t/_o-(k - q^w, «£« - iv^) 

To L ZTTVT^ V '—^ 



(A6) 



i^jT^Vv) ~ ii/s)^o(k',i£„)^t;o(k' - q,i£,„). 



k,k' 



We replace the summation over the intermediate Matsubara frequency Vs by an integration over real frequencies 
and perform the analytical continuations ie„ — > e + huj/2 + i0 and iSm — > e — fiuj /2~iQ. For the spin wave contribution 
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to the structure factor we take the hmit q, w ^ and find the expression 
72 

Im{x^CT,a(qsw,Wsw)} 
r^n^°L^(qsw;e - ?iWsw/2, w^w, 1)^ 



h tann 

2T 2r 



xRe n^"l_^(qsw;e- ?JWsw/2,Wsw,2) 



To - Tcrn^°Lcr(qsw;£ - ?iWsw/2,Wsw,0) 



E 



tanh 

27r 



xlni< X-<T,^(qsw,t^sw) 



2T 

n^°L^(qsw;e - ?ia;sw/2,a;sw,2) + 

£,Wsw)- 



r^n^'|L^(qsw;e - ?iWsw/2,Wsw, 1)^ 
To - Tcrn^"L„(qsw; £ - ?iWsw/2, i:Jsw, 0) 
To 



To - TCTn^°'_o.(qsw;£ - ?iwsw/2,a;sw,0) 



where we abbreviated as in the main text 



nl°L(q;e:w;n) 



2TTiyT„V 



E(G«(k, e + hu:/2)) (G^(k, e + ?ic./2))«(G^,(k - q, e ~ nc./2)), 



(A7) 



(A8) 



and have written the contribution by the Hikami box as 

h 



HA 



2TTiyT„V 



E(G!?(k, e)f{G^{k - q, e)) (G^,(k - q.^, s - hws^)) 



TO [2^^^ E(G?'(k, £)) (G^(k ~ q, e)) (G^,(k - q.w, e - ?.^sw)) 



(A9) 



-(^;^)'E(G^(k-s))'(G^.(k-q.w,e-7»c.,„))(G^(k',£))2(G^(k'-q,£)). 



k,k' 



The first term of (A7) is identical to Eq. (33) of the Altshuler-Aronov correction the the conductivity of a 

main text; The second term vanishes upon performing disordered metal, which is an elastic correction to the 

the integration over k. (It represents a spin-wave-induced propagation of conduction electrons and does not cause 

renormalization of the diffusion constant similar to the dephasing.) 
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